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ABSTRACT 



Edges characterize boundaries are the problem of fundamental importance in description of geometric surfaces. 
Since edge detection is in the fore front of computer visión system for detection of surfaces images as Revolution and 
Saddle surfaces needs to description, it is crucial to have a good under standing of geometric properties of Surfaces and 
edge detection algorithms. 

In this paper the comparative analysis of various surface-image edge detection techniques and the study of the 
geometric properties of Surfaces are presented. The software is developed using MATLAB 7.11. It has been shown that the 
Canny's edge detection algorithm performs better than all operators (i.e. Robert, Prewitt, Sobel, Zero-cross (Laplacian) and 
LOG) under almost all scenarios. 
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Edge detection refers to the process of identifying locating sharp discontinuities in an image like images of 
different geometric surfaces. The discontinuities are abrupt changes in pixel intensity which characterize boundaries of 
objects in an image. 

Classical methods of edge detection involve convolving the image with an operator (a 2-D filter), which is 
constructed to be sensitive to large gradients in the image while returning valúes of zero in uniform regions. There is an 
extremely large number of edge detection operators available, each designed to be sensitive to certain types of edges. 

There are many ways to perform edge detection. However, the majority of different methods may be grouped into 
two categories: Gradient and Laplacian. The visual comparison of the most commonly used Gradient, Laplacian and 
Gaussian distribution based edge detection techniques was analyzed and performed. 

In Section 2, the problem of formulation with the Gradient, Laplacian and Gaussian distribution working methods 
is presented. In Section 3, the various edge detection techniques have been studied and analyzed. Section 4 Studies the 
geometric properties of Surfaces. Section 5 discusses the advantages and disadvantages of various edge detection 
techniques by visual comparisons that have been done by developing software in MATLAB 7.11. Section 6 discusses the 
conclusions reached by analysis and visual comparison of various edge detection techniques developed using MATLAB 



There are problems of false edge detection, missing true edges, producing thin or thick lines and problems due to 
noise, etc. Suppose we have lD-image, with an edge shown by the jump in intensity (ramp edge) as shown in Figure l(a). 
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Figure 1: Image (a) with Function F(t) of Pixel Intensity Valúes 

The signal (function) for image intensity is shown in Figure l(b). The gradient (lst derivative) with respect to this 
signal is shown in Figure l(c). Clearly, Figure l(c) shows a máximum located at the center of the edge in the original 
signal. A pixel location is declared an edge location if the valué of the gradient exceeds some threshold. As mentioned 
before, edges will have higher pixel intensity valúes than those surrounding them. So once a threshold is set, you can 
compare the gradient valué to the threshold valué and detect an edge whenever the threshold is exceeded [4, 10, 16]. 
Furthermore, when the first derivative is at a máximum, the second derivative is zero. As a result, another alternative to 
find the location of an edge is to lócate the zeros in the second derivative. This method is known as the Laplacian 
(the second derivative of the signal) as shown in Figure l(d) [4, 16]. 

Now let me present the definitions of Gradient, Laplacian and Gaussian distribution in 2D-image 



GRADIENT [7, 16, 20] 

DEFINITION 2.1 In 2-dimension the gradient of function f(x,y) is: 



dx 



\Yf\=Jf?+fy , ¿-tan 1 



\fxj 



(1) 



Where I V/ I is called magnitude of gradient and 6 its direction. 

The gradient method detects the edges by looking for the máximum and minimum in the first derivative of the 
image. Let I(x, y) is the function define pixel intensity valúes of an image I (Note I(x, y) is discrete function) 



di I(x+h,y)-I(x,y) Tt , x Tt x 

— = lim^ =I(x+l ,y)-I(x, y), 

dx h^o h 



di 

— = linr 

dy h y -A) 



I{x,y+h)-I{x,y) 



(2) 



=I(x,y+l)-I(x,y) 



Where h =l,h v =1 because that / is discrete function. 
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Using pixel-coordinate notation (note: j corresponds to the x direction and i to the negative y direction, as shown 
in Figure 2) we obtain the folio wing: 
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dx 



= I(i,j+l)-I(i,j) 



Figure 2: Pixel-Coordinate Notation 
di 



=i(i-ij)-i(ij) 



(3) 



Definition 2.2 A kernel is a small matrix whose valúes are called weights. Each kernel has an origin, which is 
usually one of its positions. 

The two kernels of Gradient in direction X and y are shown in Table 1 

Table 1: The Two Kernels of Gradient 
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di di 

Now we can find — and — as follow 
dx dy 



dl_ 

dx 



= 1=1* K. 



(4) 



The gradient magnitudes are sometimes simplified by applying Manhattan distance measure as 
I V/ 1=1 I x I + I / I to reduce the computational complexity. 

LAPLACIAN [7, 16, 20] 

Definition 2.3 In 2-dimension the Laplacian of function f(x,y) is: 



v 2 / = ^4+ 



dx 2 dy 2 



(5) 



The Laplacian method searches for zero crossings in the second derivative of the image to find edges. To find the 
Laplacian of an image with pixel intensity valúes I(x, y) we first find second derivative operators oriented in the X and 

y directions. From the equations (2) we get 
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-\2t 

—i- = I(x + 2, y) - 2I(x + l,y) + I(x, y), 
dx 

|4 = I(x, y + 2)- 21 (x, y + l) + I(x,y) 
dy 2 

we can be writing (6) as 

dx 2 
d 2 I 



(6) 



dy 2 



= I(x + l,y)-2I(x,y) + I(x-l,y), 
= I(x,y + l)-2I(x,y) + I(x,y-l) 



(7) 



Using pixel-coordinate notation in Figure 2 we obtained on 
f) 1 I(i,j+l)-2I(i,j) + I(i,j-l), 



dx 2 

dy 2 



(8) 



I(i + l,j)-2I(i,j)+I(i-l,j) 



thus, the discrete Laplacian of an image with pixel intensity valúes I(x, y) is: 

V 2 I = -4I(i,j)+I(i,j+l)+I(i,j-l)+I(i+l,j)+I(i-l,j) 



(9) 



A discrete Laplacian kernel is yielded by simply adding ID second derivative operators oriented in the x and y 
directions as shown in Table 2 

Table 2: The Kernel of Laplacian by Adding the X and Y Directions 
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By also adding ID second derivative operators oriented diagonally, the most common kernel of the discrete 
Laplacian is shown in. Table 3 

Table 3: The Kernel of Laplacian by Adding the X, Y and Diagonals Directions 
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GAUSSIAN DISTRIBUTION [20] 

Because kernels are used in gradient or laplacian are approximating a first or second derivative measurement on 
the image, they are very sensitive to noise. To counter this, the image is often Gaussian filtered before applying Gradient or 
Laplacian in edge detection. 
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Definition 2.4 In 2-dimension, the Gaussian distribution function is centered on zero with standard deviation (7 
has the form: 



G(x,y)= J— e 
-42ti o 

as shown in Figure 3. 



2 2 



2cr ¿ 



(10) 




■i M 
Figure 3: Gaussian Distribution Function 

Gaussian filtering is used to remove noise from an image. The kernel of a Gaussian filter with a standard 
deviation of (7 = 1 .4 is shown in Table 4 

Table 4: The Kernel of Gaussian Filtering with (7 = 1.4 
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EDGE DETECTION METHODS 
Classical Methods 

These methods are based on gradient in edge detection. Pixel valúes at each point in the output represent the 
estimated absolute magnitude of the spatial gradient of the input image at that point. A thresholding strategy [11] is used 
for determining the local máxima in the gradient image to determine edge points. Pixels at which the magnitude of the 
gradient is above a certain threshold T can be considered to be edge pixels. The edge map E is given by 



E[m,n] ■ 



[1, \V I(m,n)\)>T 
1 0, otherwise. 



(11) 



General Algorithm [20] used in edge detection by classical operators is: 
Smooth the input image / (x, y) = I(x, y) * K c (x, y) 

I s = llx,y)*K x (x,y) 
I,=llx,y)*K y (x,y) 
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. IV/ (x,y)\=\I ¡ .\ + \I 9 \ 

• e = tan\l~II.) 

• If I V/ (x, y) \> T , then possible edge point. 

where : T is threshold , K G {x, y) is kernel of the Gaussian filtering , and K are kernels of operator edge 

detection. 

Robert's Cross Operator [12, 18] 

The Roberts Cross operator performs a simple and quick computing, 2-D spatial gradient measurement on an 
image. The operator consists of a pair of 2X2 convolution kernels as shown in Table 5 

Table 5: The Two Kernels of Robert's Operator 
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Remark: The kernel K is simply and the other is obtained by rotating K x by angle 90 . 

These kernels are designed to respond maximally to edges running at 45° to the pixel grid, one kernel for each of 
the two perpendicular orientations. The gradient magnitude is given by steps 1-4 in The previous algorithm and the angle 

of orientation of the edge (relative to the pixel grid orientation) is given by a = tan ( — ) . 

/. 4 

Prewitt Operator [7, 14] 

The operator consists of a pair of 3 x 3 convolution kernels, as shown in Table 6 

Table 6: The Two Kernels of Prewitt Operator 
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Remark: The kernel K is simply and the other is obtained by rotating K x by angle 90 . 

These kernels are designed to respond maximally to edges running vertically and horizontally relative to the pixel 
grid, one kernel for each of the two perpendicular orientations. The gradient magnitude is given by steps 1-4 in The 
previous algorithm and the angle of orientation of the edge (relative to the pixel grid orientation) is given by 

d = tan\ 1 ^-) 
h 

Sobel Operator [12, 14, 15] 

Sobel operator is similar to the Prewitt operator and is used for detecting vertical and horizontal edges in images. 
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This operator has a pair of 3 x 3 convolution kernels, as shown in Table 7 

Table 7: The Two Kernels of Sobel Operator 
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Canny Edge Detection Algorithm [5, 14] 

This operator is also based on gradient in edge detection. The Canny edge detection algorithm is known to many 
as the optimal edge detector. 

The General Algorithm [5, 20] used in the edge detection by Canny method is : 

• Smooth the image with a Gaussian filter to reduce noise and unwanted details and textures 
/ (i,j) = I(i,j)*G(i,j). 

• Compute I¡(i,j) , and then gradient of / (i, j) by using any of the gradient operators 



(roberts, Prewitt, Sobel) to get: M (i, j) = Jl¡ (i, j) + I~j (i, j) 
Threshold M: M T (i,j) 



\M(i,j), ifM(i,j)>T; 
0, otherwise 0. 



Suppress non-maxima pixels in the edges in M T obtained above to thin the edge ridges (as the edges might have 
been broadened in step 1). To do so, check to see whether each non-zero M T (Í,j) is greater than its two 
neighbors along the gradient direction 9(i, j) . If so, keep M T (i, j) unchanged, otherwise, set it to 0. 

Threshold the previous result by two different thresholds T x and T 2 (where T x <T 2 ) to obtain two binary 

images T x and T 2 . Note that compared to T x , T 2 ,they have less noise and fewer false edges but larger gaps 
between edge segments. 

Link edge segments in T 2 to form continuous edges. To do so, trace each segment in T 0 to its end and then 
search its neighbors in T x to find any edge segment in T x to bridge the gap until reaching another edge segment 



in 



2- 



LAPLACIAN OPERATOR 
Zero Crossing [7, 20] 

A thresholding strategy [11] is used for determining the zero crossings in the Laplacian image to determine edge 
points. At a zero crossing, a pixel on one side of the zero-crossing is positive and negative on the other side. The greater 
the difference between these two pixels is at a point, the more likely the zero crossing corresponds to an edge. Therefore, in 
a digital image, pixels at which this difference is above a certain threshold T can be considered to be edge pixels. 
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The edge map E is given by 

fl, L(x,y)>T 

E[m,n] = < . (12) 

[0, otnerwise. 

Laplacian of Gaussian [14, 17, 20] 

We can convolve the Gaussian smoothing filter with the Laplacian filter first of all, and then convolve this hybrid 
filter with the image to achieve the required result. Doing things this way has two advantages: Since both the Gaussian and 
the Laplacian kernels are usually much smaller than the image, this method usually requires far fewer arithmetic 
operations. The Laplacian of Gaussian (LoG) kernel can be pre-calculated in advance so only one convolution needs to be 
performed at run-time on the image. The 2-D LoG function centered on zero and with Gaussian standard deviationó has 
the form: 



2 2 

2 2 x +y 

LoG{x,y) = -\l7U(T A [\-£-^-)}e ^ 

2(7 



as shown in Figure 4. Where The x and y axes are marked in standard deviations <T . 



(13) 




Figure 4: 2-D Log Function 

The kernel of The LoG (for standard deviation <7 < 0.5) is shown in Table 8 

Table 8: The Kernel of LOG Operator with (7 < 0.5 
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Note that kernel size increases with <7 as in Gaussian filtering [8, 9, 14]. 
STUDY THE GEOMETRIC PROPERTIES OF SURFACES 

As an application, we study the geometric properties of the following surfaces 
Revolution Surface: consider a Revolution surface represented by the parametric equation 



R(u ,u ) = (u cosu ,u sinw ,e" ),u g (-3,3) 
as shown in Figure 5(a). 



(14) 
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Figure 5: Revolution Surface and its Edge by Various Edge Detection Methods 

2 

To describe the revolution surface geometric [6, 13] we define parametric lines as follows: U = const. is called 
meridians of Revolution surface. While U 1 = const. is called parallels of Revolution surface. And by calculating the 
tangents R 1 , R 2 of meridians and parallels on revolution surface we obtain 



R l = (cosw z ,smu ¿ ,e u ), R 2 = {-u l sin w 2 ,w' eos u ¿ ,0) 



(15) 



The first fundamental covariant quantities g a p and contravariant quantities g are given form g a « = < R a , 
R /3>> Sap- S fiT = SI Explicitly 



1 + e 2 " ,g n = g 2l =0,g 22 = (u l ) 2 ,g = {u l f{\ + e 2u ) 



(16) 



,11 _ §22 _ 1 „12 _ 8n _ n _ 21 22 _ £11 _ 1 

O 'O 



cr 



1\2 



8 (" ) 



(17) 



Where ^ =l + e 2 " 

The unit normal vector field A/ on (14) is given from 



{-e u cosu 2 ,-e" sinw z ,l) 



(18) 



The second partial derivatives for (14) are 



R u = (0,0, e" ), R l2 = (-smu z ,cosu ¿ ,0), R 22 = (-u l cosw 2 ,-w' sinw 2 ,0) 



(19) 



The second fundamental quantities a P ' take the form 



7 '^12 — ^21 — 0,-Zj 



1 u 1 

u e 



22 



(20) 



Gaussian curvature K and mean curvature H and principal curvatures k 1 , k 2 are given from 
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L _ L n L 22 L 12 



2u l 



K - A, A, - - . 



_ 2 1 fl 

S <?11<?22 #12 U fel 



(21) 



2H = k l +k 2 = g a/í L a g = g l % 1 +2g u L l2 + g 2 % 2 = 



uV 2 



(22) 



Where 7] l = 1 + U l + e 2 "' . From the equations(21), (22) we get 



..i r 



K = 



2Í 



2] 



i/2 



(23) 



Where £ = +é> 2 " 

Saddle Surface (Ruled Surface): consider a saddle surface represented by the parametric equation (as shown in 
Figure 6(a)) 



R(u\u 2 ) = (-O 1 +u 2 ),-(u 2 -u l ),u l u 2 ) 
2 2 



(24) 



To describe the Saddle surface geometric [6, 13] we define parametric lines as follows: U = const. is called 
meridians of Saddle surface. While w' = const. is called parallels of Saddle surface. And by calculating the tangents 
R l , R 2 of meridians and parallels on Saddle surface we obtain 



„ A 1 2x n ,11 K 

^ =(---,« ),R 2 = (-,-,u ) 



(25) 



The first fundamental covariant quantities g a p and contravariant quantities g are given form g ( 
>. g a p-8 Pr = S y a Explicitly 

= ^ + (" 2 ) 2 '<?i2 = £21 =u l u 2 ,g 22 =^ + (u l ) 2 ,g = -£ 2 



= <R„ 



(26) 



where £ 2 = 2(u 1 ) 2 + 2(u 2 ) 2 + 1 
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Figure 6: Saddle Surface and its Edge by Various Edge Detection Methods 

u = 2+4CK 1 ) 2 12 = -4aV = 21 22 = 2 + 4( M 2 ) 2 

b 2 b 2 b 2 

The unit normal vector field N on the saddle surface (24) is given from 

N(u\u 2 )= <Rl, l ^ 2 > =-¡={-u l -u\u 2 -u\\) 



The second partial derivatives for (24) are 

R n = (0,0,0), R 12 = (0,0,1) = R 21 , R 22 = (0,0,0) 

The second fundamental quantities L aB =< R aB , N > take the form 



Ljj — 0,L 12 — L 2l — j— ,L 22 — 0,L — 



Gaussian curvature K and mean curvature and principal curvatures k x , k 2 are given from 

L -4 

K = Lk 7 = - = — *0 

12 g a 



2H=k ] +k 1 = g ap L„ R =- 



O 1 2 



3/2 
2 



by solving the two equations(31), (32) we get 

4w^ 

Ti 



A 1 2 

k - l\t - k --2ÍF 

"1 — z "Vb2 >-3í2 ' 2 — Z "V'2 



4mV 

Í-3Í2 
b2 



where £ = , 7? 2 = (1+ 2( M 1 ) 2 )(1+ 2( M 2 ) 2 ) 
b 2 



(27) 



(28) 



(29) 



(30) 



(31) 



(32) 



(33) 
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ADVANTAGES AND DISADVANTAGES OF VARIOUS EDGE DETECTION TECHNIQUES 

Edge detection of all two types (Revolution and Saddle surfaces) was performed on the above Figures 5(a)-5(g) 
and 6(a)-6(g). Canny yielded the best results. This was expected as Canny edge detection accounts for regions in an image. 
Canny yields thin lines for its edges by using non-maximal suppression. Canny also utilizes hysteresis with thresholding. 

As edge detection is a fundamental step in computer visión, it is necessary to point out the true edges to get the 
best results from the matching process. That is why it is important to choose edge detectors that fit best to the application. 
In this respect, we present some advantages and disadvantages of edge detection techniques [14, 16, 19] as follow 

• Classical Methods (Roberts, Prewitt, Sobel) 

• Advantages: Simplicity, Detection of edges and their orientations. 

• Disadvantages: Sensitivity to noise, Inaccurate. 

• Zero-Crossing (Laplacian) 

• Advantages: Detection of edges and their orientations. Having fixed characteristics in all directions. 

• Disadvantages: Responding to some of the existing edges, Sensitivity to noise. 

• Laplacian of Gaussian(LoG) 

• Advantages: Finding the correct places of edges, Testing wider área around the pixel. 

• Disadvantages: Malfunctioning at the corners, curves and where the gray level intensity function varíes. 
Not finding the orientation of edge because of using the Laplacian filter. 

• Gaussian (Canny) 

• Advantages: Using probability for finding error rate, Localization and response. Improving signal to noise 
ratio, Better detection specially in noise conditions. 

• Disadvantages: Complex Computations, False zero crossing, Time consuming. 
CONCLUSIONS 

Since edge detection is the initial step in object recognition (as surfaces geometric), it is important to know the 
differences between edge detection techniques. In this paper we studied the most commonly used edge detection 
techniques of Gradient-based and Laplacian-based edge detection for geometric surfaces (Revolution surface and Ruled 
surface) as shown in Figures 5(a)-5(g) and 6(a)-6(g). The software is developed using MATLAB 7.11. 

Gradient-based algorithms such as the Prewitt filter have a major drawback of being very sensitive to noise. 
The size of the kernel filter and coefficients are fixed and cannot be adapted to a given image. An adaptive edge-detection 
algorithm is necessary to provide a robust solution that is adaptable to the varying noise levéis of these images to help 
distinguish valid image contents from visual artifacts introduced by noise. The performance of the Canny algorithm 
depends heavily on the adjustable parameters, which are the standard deviation for the Gaussian filter, and the threshold 

valúes, T x and T 2 . also controls the size of the Gaussian filter. The bigger the valué, the larger the size of the Gaussian 

filter becomes. This implies more blurring, necessary for noisy images, as well as detecting larger edges. As expected, 
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however, the larger the scale of the Gaussian, the less accurate is the localization of the edge. Smaller valúes ofó imply a 
smaller Gaussian filter which limits the amount of blurring, maintaining finer edges in the image. The user can tailor the 
algorithm by adjusting these parameters to adapt to different environments. Canny's edge detection algorithm is 
computationally more expensive compared to Sobel, Prewitt and Robert's operator. However, the Canny's edge detection 
algorithm performs better than all these operators under almost all scenarios. Evaluation of the surface images showed that 
under noisy conditions, Canny, LoG, Laplacian, Sobel, Prewitt, Roberts's exhibit better performance, respectively. 

Nassar et.al. studied the image processing problem using the ridge and Ravens on the surface through the 
máximum and minimum of the principal curvature [1, 2, 3]. In a forthcoming paper we will find method to edge detection 
of surface-image based on geometric properties of surfaces, where we will try to overeóme flaws in previous methods. 
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